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Abstract 

We consider quantum interferences of ballistic electrons propagating inside fractal structures 
with nanometric size of their arms. We use a scaling argument to calculate the density of states of 
free electrons confined in a simple model fractal. We show how the fractal dimension governs the 
density of states and optical properties of fractal structures in the RF-IR region. We discuss the 
effect of disorder on the density of states along with the possibility of experimental observation. 

PACS numbers: 05.60.Gg Quantum transport. 73. 23. Ad Ballistic transport 78.67.Bf Optical properties of 
low dimensional, mesoscopic, and nanoscale materials and structures. 
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I. PECULIARITY OF METALLIC NANOFRACTALS 



Ramified structures are widely observed in nature at scales from the microscopic world 
up to the human size. They have been studied in various contexts and in different domains 
of science: biology, physics, chemistry, etc. Surface Science is one particular field where the 
ramified semi-metal 1 ^, semiconductor, metallic^ or dielectric^ structures may range from 
the nanometric up to the micrometric sizes. The mean free path of electrons in metals is 
usually of the order of 10 2 -10 3 nm depending on the kinetic energy. Therefore electrons prop- 
agating ballistically in metallic nanostructures may manifest essentially quantum behavior 
associated with strong interference of their De Broglie waves in contrast to the diffusive^ or 
hopping 9 behavior intensively studied during the last decades. The combination of quantum 
ballistic motion and ramified geometry suggests to consider the interference of electrons in 
a fractal metallic structure confining their propagation. 

Tree-like structures is a natural example of fractals. Results obtained for quantum par- 
ticles moving on tree-like latices^, for the quantum localization in the framework of sparse 
random matrix models^ topologically similar to trees, and for quantum systems with tree- 
like hierarchy of interactions^ have revealed a certain universality associated with such a 
topology, that persists in different physical situations. Therefore for tree-like fractals one 
can also expect a universality of the quantum properties related to their specific geome- 
try. Moreover, the key property of fractal structures is the invariance under certain scaling 
transformations. Therefore considering quantum dynamics of electrons on fractal trees we 
take advantage of the scaling arguments 13 . Note that it is equally important to study the 
properties of ensembles of isolated or interacting fractals placed together at a surface, since 
it is experimentally difficult to address a single nanometric object. Models of such ensem- 
bles might be also of interest for consideration of conductivity of thin films^, heterogeneous 
catalysis of nanometer larger silver particles^, quantum dot networks^, and in other do- 
mains. 

In this paper, we consider the simplest tree-like fractal with identical length of the 
branches at each generation and symmetric nodes as a support of ballistically propagat- 
ing electrons. We introduce a single geometrical parameter a which gives the ratio of branch 
lengths for successive generations. We shall see that this parameter is closely related to the 
fractal dimension of the tree. We show that the density of the one-electron states manifests 
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a power law dependence on the momentum near zero energy with the power index being 
the fractal dimension. It is consistent with the result^ for the low momentum asymptotic 
of Green functions in systems of fractal dimensionality. Note that this property is typical 
of fractals since linear objects of the same size do not have quantum states close to zero en- 
ergy according to the Born-Sommerfeld quantization rule. We demonstrate the macroscopic 
manifestations of this power law in optical properties of surfaces covered by the nanometric 
ramified structures by calculating the reflectivity in the RF-IR frequency domain. Finally, 
with the help of a simple random matrix approach^ we consider the role of irregularities 
in fractal structures resulting from the statistical distribution of branch lengths and nodes 
asymmetries, that does not require to allow for the level-level correlations in the ballistic 
regime. 

We formulate the problem in terms of the Green functions of a particle propagating 
along the fractal. We employ the momentum variable which is natural for consideration of 
the interference phenomena, whereas the energy dependence is given by the dispersion law 
E = E{jp) specific for each type of systems. It allows one to implement the results for any 
particular dependence of the particle energy on the momentum which are usually different 
for metals and for semiconductors: for a free particle E = p 2 /2m, where m is the mass of the 
particle, whereas for metals E = Vf\p\, where Vf is the Fermi velocity. One-particle Green 
functions are obtained following the standard quantum field formalism widely developed in 
various textbooks 19 . Quantum state density g(p) and several other properties such as linear 
dipole response R(u) or conductivity <j{uj) at a frequency uj can be found with the help of 
its retarded G R (E) and advanced G A (E) Green operators via the relations 20 



where d is the dipole moment operator, j is the current operator, and p is the density matrix. 
By the analogy to a photon propagating in a Fabri- Perot resonator, we can take into account 
only the coordinate parts G(p, x, x') of the Greens operators at a given energy E(p) ignoring 
the resonant denominators (E — E(p) + zO). The latter can be factored out during the 
consideration of the interference phenomena and have to be restored only at the last stage, 
prior to substitution to Eqs.tJTJ). Note that in the case of ballistic propagation the coordinate 



g(p) = - 

R(u) = Tr G A \ 
a(u) = Tr G A 



-±Im Tr G R {E) 
(E)dG R {E + u)dp{p), 
(E)3G R (E + u)dp(p) 



(1) 
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part of the product G A (E)G R (E + uj) can be written in a single factor G(p) depending only 
on the momentum p = u>.p' E , where p' E = associated with the energy shift uj. For g(p), 
Im R(ui), Re <j(uj) the allowance for denominators yields the Dirac 5-functions of energies 
which disappears after taking the trace. Therefore these parameters responsible for the 
absorption of electro-magnetic radiation can be calculated directly when we replace G G R 
in Eqs.Q by G(p). The Kramers-Kronig relation then yields the dispersive parts Re R(uj) 
and Im a(u). In this paper we therefore call "Green function" the coordinate part G(p, x, x') 
of G A ( Pf )G R (p + p f ). 

For metals the density matrix is given by the Fermi step f)(p) = Vfn e Q(—p) where n e is 
the electron state density in metal near the Fermi surface and the Fermi momentum is taken 
as a reference point. The dipole moment operator d in the momentum representation reads 
d = ied/dp where e is the electron charge (we set h — 1), whereas the current operator j is 
simply pe/m. Therefore Eq.Q takes the form 

g(p) = i Tr G(p) 
Im R(u) = v f e 2 n e || =o Tr G (p + ^) (2) 



r.- 



Re a(u) = Tr G , 



where we have taken into account the relation dQ(p)/dp = —5(p). The trace operation now 
implies only summation over all closed trajectories in the coordinate space corresponding to 
a given momentum in complete analogy with the Fabri-Perot resonator. 



II. THE MODEL OF FRACTAL 

We model a fractal by three trees with trunks joint in a node at the fractal center (FigJTJ). 
Each of the trees starts with a trunk of length L and is built by recursive attaching at each 
terminations two homothetical branches scaled by a factor a. The homothetical factor a is 
the main parameter of the model. It governs all geometrical properties and in particular the 
fractal dimension which is the main physical parameter. For a > 1 branches are longer at 
each step, whereas for a < 1, branches are smaller as n increases, which is always the case 
in our consideration as we shall see. Electrons propagate ballistically along the trunks and 
branches until they reach a node where three branches are attached symmetrically at the 
angle 2tt/3 as shown in FigO Nodes scatter the electrons backwards and forward into the 
attached branches. 
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A. Nodes model 



The branches joining a node have different length which depends on the index n numerat- 
ing the generation, that is the number of nodes which separates the branch from the fractal 
center. Two branches are of the length L n = a n L whereas the branch closest to the trunk 
has the length L n -\ = a n ~ 1 L. If we stop the development of the tree at a given n = N, the 
last rightmost branches have a length = a N L and the total number of such branches is 



FIG. 1: Model tree consisting of a trunk of length Lq = L and successive adjunction of branches. At 
each generation n, two branches of length L n = a.L n -\ are attached to the previous branch. The 
scaling parameter a governs the tree morphological properties, mass, length and fractal dimension. 
We restrict to 1/2 < a < 1. The scattering matrix (Eq0J) couples incoming fluxes to outgoing ones. 
The fractal of radius Z is built by attaching three identical trees to the root. 

Having arrived at a node an electron either scatters into the two attached branches 
with equal (due to the symmetry) probability or returns back with a different probability. 
The node is formally described by a unitary 3x3 scattering matrix S with the matrix 
elements Sjj/ coupling three outgoing probability amplitudes 4>j of the electron to the three 
incoming ijjj ones, where the marker j assumes the values, I, r, and b for the left-scattered, 
right-scattered, and the back-scattered amplitudes, respectively. The relation among the 
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amplitudes reads 



\4nj 



' S bb Sbr S U \ ( lfj b » 



V 



(3) 



S r b S rr S r [ 
Sib Si r 8 u 

Apart of the unitarity, the matrix S should satisfy two more requirements imposed by 
the node symmetry and by the long wave limit. The symmetry requirement implies that 
the probability amplitudes for the left- and the right-scattering given by the coefficients s r b 
and sib respectively, are equal. Moreover, the symmetry with respect of the node rotation at 
the angle 27r/3 implies that all other off-diagonal coefficients also have the same value. We 
also assume that no quantum defect is associated with the scattering at the node. In the 
long wave limit it implies that no phase shift is introduced during the scattering process, 
and hence all the parameters Sjji are real. These three requirements together yield 
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(4) 



as the only choice for the scattering matrb4=i 



B. Scaling factor and the fractal dimension 

Now we relate the typical length L of the system and the scaling factor a with the fractal 
dimension employing the self similarity aspect of the problem. In fact, in the general case L 
is not the only typical length scale in the problem. The homothetical factor a governs most 
of the advanced morphological properties of the model tree. The total length of the tree 
with truncated branches of iV + 1-th generation reads 

N N 1 _ n N+l 

z N = Y,L k = J2 akL = \_ a ( 5 ) 

fc=0 fc=0 

This expression imposes a first limit on a : for < a < 1 the length Z^ converges to a finite 
value 

Z = ^~ (6) 
I — a 
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whereas for a >z 1 it diverges. We consider the fractals of a finite size only. Actually, the 
radius of a tree is given by a more complicated expression and should take into account 
the geometrical arrangement of the branches with 27r/3 angle between them. The exact 
calculation for the diameter gives = L(2 + o)(l — a N+1 )/(l — a 2 ) which also converges 
when N — > +00 for a < 1 to the value D = L(2 + a)/(l — a 2 ). 

The mass Mjv of the tree that is the sum of the lengths of all branches is given as 

N N 1 (?n) N + 1 

M N = T 2 k L k = V 2 k a k L = (2fl) L (7) 
/ — ' ' 1 — 2a 

k=0 k=0 

which converges to M = L/(l — 2a) for a < 1/2 and diverge for a > 1/2. We are interested 
in the regime where the mass of the fractal is infinite, and therefore a ranges from 1/2 to 1. 

The model fractal has the same fractal dimension as its consisting trees. The fractal 
dimension of a tree is given by a standard evaluation 22 ' which is now widely used. It implies 
the calculation of the minimum number N(e) of disks of diameter e needed to completely 
cover the whole tree. In a fractal structure, gradual decreasing of e reveals new details caus- 
ing N(e) to vary non trivially as e~ Dh where Dh defines the so called Hausdorff-Besicovitch 
fractal dimension. 

Let us implement this definition in our case of tree-like fractal. In order to find the number 
N(e) of e-sized disks required for covering the tree we make use of the scaling arguments. 
Let us take the infinite tree and applying to it the homothetical factor a. One obtain another 
tree which also has the same infinite structure but starts with a smaller trunk of length aL. 
This a-contracted tree can be considered as an element of the original tree, namely its first 
generation branch with all the branches of subsequent generations attached. The size of the 
discs covering this branch is apparently a times smaller compared to original discs of the 
radius e. When we attach two a-contracted trees to a trunk of length L we recover the 
original form of our fractal with branches covered by 2N(e) discs of radius ae. One requires 
L/ae additional discs to cover the trunk. We therefore obtain the equation 

N{ea) = 2N(e) + — (8) 

CIS 

determining an asymptotic behavior of N(e) for e — > 0. 

We look for the solution of Eq.fjHJ) in the power-law form N(e) ~ e~ a with a > 1. It 
implies that the second term in the right hand side of Eq.flHJ) can be omitted with respect 
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to the first term as e — > 0, and we arrive at (ae) a = 2(e) a + o(e a ). It yields 



a 



In 2 
In a 



(9) 



which is the well-known Hausdorff-Besicovitch fractal dimension of a self similar recursively 
built fractal^. Equation © gives fractal dimension greater than 1 in the case a > 1/2 
corresponding to an infinite mass M. We also restrict ourselves to the case a < 1 corre- 
sponding to a finite size Z of fractals. In this regime the spectral peculiarities typical of 
such structures manifest themselves in the most interesting way. 
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FIG. 2: Fractal dimension as function of the scaling parameter a. For < a < 1/2 the tree has 
a finite mass and its fractal dimension is dominated by the trunk : = 1. For 1/2 < a < 1 the 
mass is infinite whereas its length remains finite : = a = — ln(2)/ln(a) 



III. GREEN FUNCTIONS AND QUANTIZATION OF THE FRACTAL STATES 

The Green functions G(p) generally given by the Feynmann path integral can be found 
for the particular case of a tree-like fractal structure from recurrent relations formulated in 
terms of the Green function of a free one dimension particle G (p) = Q(x—x') exp[ip(x — x')] 
propagating along the branches and the scattering conditions Eq. (jHJ) at the nodes. We derive 
a recurrent relation for these functions and thereby determine the spectrum of the eigen state 
density. 
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A. Recurrent relations 



The idea of derivation of the recurrent relations is illustrated in Fig|3^). By X n (p) 
we denote the unknown exact Green function for the particle leaving a chosen node of ra- 
th generation and returning back after multiple scattering in all the variety of nodes of 
subsequent generations connected to the chosen node. Then the Green function of 
the previous generation can be considered as a result of the free propagation of the particle 
towards the ra-th node followed by the multiple scattering at this node resulting in the direct 
back scattering and the scattering to the attached branches followed by the multiple 
returns and back-scattering in the nodes and branches of the subsequent generations. One 
finds the result of all these multiple scattering events by considering the relation Eq.fjSJ) 
among the incoming <fi and outgoing ip amplitudes with the allowance for the fact that they 
are related by the condition 

VV = X n{p)<f>i,r (10) 

which holds by the definition of Green functions. 
The free propagator G a (p) gives the relation 

•0 = exp[ipLa n ](p b ; ip™ = exp[ipLa n ]4> (11) 

between the amplitudes ipb and <pb of the waves incoming to and outgoing from the node 
n along the branch attached to the node n — 1 and the amplitudes <fi and ip of the waves 
outgoing from and incoming to the latter. Here we do not specify whether ip m > out corresponds 
to the right scattered or to the left scattered amplitudes at the node n — 1 since the relation 
are identical for both cases. The scattering matrix Eq.(0J) and the condition ip = X n _i<f) 
together with Eqs. (j3llOHH) yield the exact recurrence relation for the Green functions 

*n-i(p) = exp[2^La"] 1 ~ 3 f" (p) (12) 

X n {p) - 3 

Equation ()12j) maps the Green function X n (p) of a n-th node to the Green function 
X n _i{j>) corresponding to a node of the previous generation. As we are interested in the 
high n behaviour, this equation has to be inverted to obtain the expression of X n (p) as 
function of X n _i(p). Changing the n index to n + 1 we have 

exp[2^La" +1 ]+3X n 
An+l[P) ~ 3exp[2zpLa«+i] + X n (p) {L6) 
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FIG. 3: Recurrent relations for the Green functions, (a) A Green function Xn_i attached to a 
parent node consists of a free propagation Go followed by the scattering in the node of the next 
generation and back propagation. Amplitudes by the outgoing and incoming waves in each branch 
differ by the factor Xjv-i- (b) The mapping corresponding to Eq. (|l«-Sj) in the long wave limit p = 
has two stationary points. One corresponds to a regular back-scattering with the phase shift — tt, 
whereas the other does not yield any phase shift and gives rise to the essentially fractal domain of 
the spectrum near zero energy, (c) Boundary conditions at nodes corresponding to the stationary 
point —1 (left), and +1 (right). 

that we make use in FigJ3)D) for p = 0. 

This mapping EqJ13l has two stationary points X st = =pl. Both of them have physical 
meaning. The negative sign corresponds to the regular situation when the reflection of 
the wave function from a node occurs with a phase shift — ir exactly in the same way as 
the reflection from an infinite vertical barrier implied by the boundary condition ip = 0. 
The positive sign corresponds to a free border ip' x = when the wave goes through the 
node and returns back with no phase shift, as shown in FigOb). The latter case changes 
the quantization rule for a particle moving in a branch confined by such nodes from both 
sides allowing the eigen states at zero energy that do not exist for the regular confinement. 
Vicinity of this stationary point gives rise to a specifically fractal domain of the energy 
spectrum at small values of the energies and momenta. 
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B. Scaling 



Now we make use of the scaling arguments and find the Green function in the long wave 
asymptotic and large n. The scaling assumption implies that X n _i(p) = X n (ap), which 
means that the Green functions X(x) corresponding to the branches of any generation are 
functionally identical and differ only by scaling of the argument x = a n p. Therefore in the 
long wave asymptotic where exp(2IpLa n ) —>■ 1 Eq. (jl2J) takes the form 

of a functional equation, where we have employed a small dimensionless argument x = Lp 
instead of p. 

This equation has an exact solution 

1 — r a 

X(x) = (15) 

with a given by Eq.® and yields an asymptotic expression 

_ , * 1 - (a n Lp) a 

XJLp) = 4-. 16 

Equation f)16|) holds for small arguments. However, even for a large values of x = Lp or 
small n an accurate numerical approximation can be obtained with the help of few iterations 
of the exact recurrent relation Eqll2l For low p and for a = 2/3, with EqE3as a starting 
point, say n = 10 iterations of Eq|T21 gives a good approximation within 1% compared to 
the exact solution EqfTBl 



C. Quantization and state density 

Now we are in the position to perform the quantization of the particle motion on the 
entire fractal and determine the density of the energy eigen states. For the purpose we 
consider the root node at the center of the fractal with three trunks attached and calculate 
contributions of all closed trajectories that start and end in a point of one of these trunks 
close to the node. The trajectory sum starts with the zero length trajectory which gives 
the contribution 1. The trajectory first going to the trunk and returning back gives the 
contribution X (Lp), whereas the contribution of the trajectory which first goes to the node 
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is B = (1 — 3X )/ (Xq — 3) according to Eq.([T2|). The trajectories of the second order give 
XqB and BXq whereas the third order results in XqBXq and BXqB. The overall sum reads 

Tr G(p) = 1+X + B + X B + BX + .. 

_ (1+X )(1+B) = 2 1+X Q ' ' 

1-BXq 3 1-Xq ' 

as it follows from summation of the geometric series. 




FIG. 4: Density of states for different fractal dimensions a = 1.18 corresponding to a = 5/9 (solid 
line), a = 1.41 corresponding to a = 11/18 (dash-dot line), and a = 1.71 corresponding to a = 6/9 
(dashed line), calculated with the help of Ea. (|17|) where Xq has been obtained from Ea. l|16jl after 
10 iterations of Ea. H12j) . The fractal diameter is the same for all fractal dimensions. Vertical lines 
shows positions of the levels in a one dimensional potential well of a width equal to the fractal 
diameter. 

In the long wave limit, injecting Eq. (jl5|) into Eq. (jl7|) we find 

Tr G(p) = l(Lp)- a , (18) 

which shows that at small energies the density of fractal energy eigen states follows the power 
law dependence on the momentum with the power index given by the Hausdorff-Besicovitch 
fractal dimension. In Fig0]we illustrate the difference between the fractal spectrum found 
from Eq. (jl7|) and the spectrum of a one dimensional particle moving in the potential well 
of the width 2Z — 2L/(1 — a) suggested by Eq.© for the fractal diameter. One clearly 
sees that the fractal boundary conditions at the nodes corresponding to the stationary point 
X = 1 of mapping Eq. (j!4)l result in the appearance of the spectrum near zero energy, where 
the potential well does not have eigenstates. 
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IV. NANOFRACTAL RESPONSE TO IR-RF FIELD 



Let us consider now the optical response of the nanofractals calculating the reflectivity of 
a transparent support surface covered by fractals as a function of the incident field frequency. 
We start with the case of isolated fractals each of which independently contribute to the 
reflectivity. The typical frequency domain can be estimated as the inverse of the typical time 
of flight of an electron across the fractal given by the Fermi velocity divided by the fractal 
size Z Eq.(J£)J), which for the fractals^ of lOOnm corresponds to the THz frequency domain 
that is far IR or short RF radiation. Then we consider the case of "merging" fractals, when 
the neighboring fractals irregularly placed at the surface can interact with capacitor-like 
connections via their most closely approaching terminations. 

A. Isolated nanofractals 

The Maxwell equation 



for a plane electromagnetic wave E incident normally to a surface covered by isolated frac- 
tals at x = allows one to find an intensity of the reflected field E r provided the specific 
conductivity <t s (uj) and the specific dipole susceptibility R a (u) of a unit surface area are 
known. The inhomogeneities of the surface have to be much smaller compared to the wave- 
length of the wave and the thickness b of the fractal layer. For a wave incident at an angle 
to the surface the same equation is valid for the tangent component of the field, whereas the 
normal component is not affected by the layer of the fractals. The continuity condition for 
the tangent field and the jump of its derivative across the surface 




— <*(-) iua s (tu) + — R s {u) E 




(19) 



a (E - E P - E t ) 



E + E r — Ej 



(20) 




E c 2 +2ir[ico a {u))+u)R s (u))\ 



(21) 



~ =^[ica s (cu)+ujR s (uj)} 
for the ratio of the reflected and the incident field amplitudes. 
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Equations (|2TT5j) yield 



Im R(u) = -N f avj^K^)- 

2e 2 n e uj at I Lu> \ —a 



Re aJw) = ^^N f v 



(22) 



where iVy is the number of fractals per unit area. We replace the product Nfn e by the specific 
density of states n s of the fractal material near the Fermi surface multiplied by the total 
volume V of the material deposited per unit surface, express the trunk size L = Z (1 — a) in 
terms of the typical fractal size Z Eg. ([5)1 and the fractal dimension a Eq.(jOJ), and substitute 
a/vfX instead of the product n s e 2 where a is the residual conductivity^ and A is the electron 
mean free path in bulk metal. In the last replacement we assume that n s = N e i/pfVf where 
N e i is the density of the metal electrons. We arrive at 

Im RJu) = ^EZ_(^)-« 

Non-analytical behavior of these dependencies at uo = does not allow one to determine 
the dispersive parts Re R s (u) and Im cr s (u;) from the Kramers- Kronig relations. However 
the latter should be of minor importance provided the transparent material supporting the 
fractal at its surface has a refraction index r different from 1. In the latter case 



E r VbaSm{Zu/v f )- 



\ a 



avf ijj 



Pf v f 



1 — r 

+ z, • (24) 

1 + r v ' 



E cA 3 (2 1 /" - 1) 

The simplest possible way to find the missing parts is to take an analytical continuation of 
Eq. (|2*3*j) to the complex plane such that giZiojvf) ~ co>~" vanishes at the negative part of 
the real axis. 



B. Ensemble of nanofractals 

When the size of the fractals becomes larger than the inter-fractals distance, the model of 
isolated fractals fails, since the dipole approximation for the response is not any longer valid. 
In the same time, allowing for the contribution related to the conductivity we have to take 
into account the points of the closest approach of neighboring fractals, where the potential 
difference experience large changes. These domains work as capacitors that assume the 
main part of the dipole activity of the system. When the ramified structures are randomly 
distributed on the surface but not yet result in the electric current percolation, as it is the case 
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for the experimental work^ for instance, the fractal ensembles conform the Dykhne model^. 
Formulated for a two-phase random conducting surface with the conductivities o\ and 02 
different for different phases this model yields the macroscopic conductivity a e ff = ^Jo\02 
which immediately suggests 

1/2 



eff 



AaVuo 



3Xmvj (2 1 /" - l) c 



t iub 



(25) 



for the effective conductivity of the fractals covering the surface. Here we have assumed 
that the capacitors of a plate size b separated by a mean shortest inter-fractal distance d are 
subjected to the potential difference ZE accumulated on the distance of the fractal size Z. 
Substitution of Eq.^ to Eq.flHJ) yields 



E r 
E~ 



.2irub 2 



v f J 



1 



(26) 



cv f l6XmZd (2V a - 1) Q 

where we have omitted the real part of the effective conductivity as small relative to the 
support contribution. 



C. Random fractals 



Thus far we have been considering the model of an ideal fractal with a high symmetry 
and an exponential variation of the branch lengths with generation number. In order to get 
an idea of how close can be such a model to the reality we now consider an ensemble of 
irregularly distorted fractals. The simplest way to model the random distortion is to treat it 
as a perturbation of the fractal Hamiltonian by a random matrix with a given mean square 
{V 2 ) of the matrix elements. The transformation rule 

G(E) = G(E(E)) (27) 



E = E + (V 2 ) TtG(E) (28) 

suggested by one of the authors 18 as a simple way to solve the Pastur— equation describing 
such a perturbation. Eq. (|2Tj) relates the ensemble averaged perturbed Green function G(E) 
with the unperturbed one G(E) depending on a transformed argument E(E). The trans- 
formation E(E) follows from the solution of a nonlinear algebraic equation (Eq. ()28j) ) which 
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allows one to find E for each E selecting from many possible solutions the one continuously 
changing from — oo to oo for E varying in this interval. By the same replacement of the 
argument one can obtain all other linear properties of the randomly perturbed system. 

Comparing Eqs. (jl|2j) and Eq. ljlH)) with the allowance of the condition g(E < 0) = one 
finds an expression 

n e<£ ) = ^mRF^^ < 29) 

consistent with the state density Eq.Q for both the positive and the negative energies. 
The constant go enters as a cofactor of the other unknown quantity (V 2 ) and both factors 
together form a single energy parameter W = go(V 2 ) responsible for the strength of the 
random perturbation. We substitute Eq. (jHEEHI) to Eq. (}28|) and obtain 

-Z(l - 2- l ' a )E 



n ~ W 

E = E + 

sni7ra 



(30) 



One sees that by introducing an energy scaling factor E = 8F with F = 
(Z(l - 2- 1 / a )/v f ) a/{a+1) ty-V(«+i) equation flflfl can be reduced to the form 

8 = 8 + (31) 

Sin7TQ! 

which does not contain parameters other than the fractal dimension. 

In order to find the universal dependencies 8 (8, a) there is no need to solve the equation 
(131 j) . After the replacement 8 = —Ke l9 , one eliminates k employing the fact that 8 is real 
and finds the dependence 8(8) in a parametric form 

e(e) = -e i0 

v ' \ sm Ssin7ro / (32) 

8(6) = -( r^ nad )^ cos e+ ( r s j nad cosad 

V / V sm osmira I \ sin osmna I 



sin7ra 



The imaginary part of 8(8, a) /it shown in FigJSJas a function of the energy 8 for different 
fractal dimensions a yields the shape of the state density g a ,w{E) = Flm8(E/F, a) /it which 
for the case of irregular fractals should replace the factor TrG(p)/7i = -^[Zcu/vf(2 1 ^ a — 
1)]~ Q ~ 8~ a in the expression Eq. (|^jl as well as in Eqs. (j23j) for the dipole response and the 
conductivity and in Eq. (|2f)j) for the effective conductivity of a disordered surface. It yields 



Ej v I cA 



Im- 

E 2-Kuib 2 



ELA. j L_ 

c du> p f v f 



cv 



f 



XmZd 



( Zuj/vf 
9a,W I 2 l/a_i 



(merging fractals) 



(33) 



for the absorption of isolated and merging fractals. 
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FIG. 5: Universal forms of the quantum state density profiles for randomly perturbed fractals with 
the fractal dimensions a = 1.18 (solid line), a = 1.41 (dash-dot line), and a = 1.71 (dashed line). 
These profiles do not depend on the size of the perturbation which results only in scaling of the 
energies. 

V. POSSIBILITY OF OBSERVATION 

We conclude by discussing the possibility to observe the optical manifestations typical 
of fractal structures experimentally, for realistic parameters of nanostructures. We take 
PfVf ~ 5eV, Vf/c ~ 1CT 2 for the Fermi velocity and momentum, V ~ Inm for the mean 
thickness of the fractal material at the surface, b ~ lnm for the cross section size of the 
fractal branches, Z ~ A ~ lOOnm for the fractal radius of the order of the mean free path 
on an electron in metal, a[Ag]/sQ = 6.3 10 7 /8.85 10 -12 sec _1 for the silver bulk conductivity 
in CGS units, and d ~ lOnm for the inter fractals distance. For the frequency u [THz] we 
take the units 10 l2 Hz natural for the electrons moving inside the nanometric sized objects. 
In order to be specific we chose the fractal dimension a = 1.41 which corresponds to the 
scaling factor a = 11/18. In this regime from Eq.(|33|) one finds 



which corresponds to the energy absorption at the level of 10 -4 . Such a small absorption is 
associated however with a phase shift of a few degrees, which is normally detectable by the 
ellipsometric measurements in the optical domain. The same estimate also can serve as the 
detection limit for IR domain whereas the internal reflection technique should be even more 
sensitive. 




-10 



2 



+ ^ 10 ~ 2 ] 9a,w (0.71 u) isolated 
5^ [ g a ,w (0-71 oj)} 1 ^ 2 merging 



(34) 



17 




FIG. 6: Density of states (a), optical response of surfaces covered by isolated fractals (b) and, 
merging fractals (c) when the inter-fractal distance is smaller compared to the fractal size. For 
a regular fractal (solid lines) one sees the power law dependencies corresponding to the fractal 
dimension a = 1.41 chosen. No typical energy reference exists for the unperturbed fractals, whereas 
for disordered fractals the typical energy is given by the disorder parameter W = go < V 2 > which 
is small (dashed line) medium (dotted line) or large (dash-dot line) with respect to the energy unit 
chosen. 

The dependencies Eq. ()34)1 are shown in FigElfor different sizes of the disorder parameter 
in the regime of both isolated and merging fractals. The power law dependence correspond- 
ing to the ideally symmetric fractals manifests itself as an asymptotic dependence for the 
irregularly perturbed fractals when the frequency exceeds the typical size of the parameter 
W governing the disorder. 
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